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ABSTRACT 

A review of the status of lattice simulations in particle physics phenomenology is presented. Recent computations of 
leptonic decay constants of light and heavy mesons, and of the Isgur-Wise function relevant for semi-leptonic decays 
of _B-mesons, are discussed in some detail. Calculations of other quantities are briefly outKned. The systematic errors 
inherent in lattice simulations, and procedures to reduce and control them, are described. 



1 Introduction 

In this talk I will briefly review some of the important 
questions in particle physics phenomenology which are 
being studied using the lattice formulation of (Quantum 
^T^Chromodynamics) QCD and numerical simulations. I 
ONhope to convince you that "Lattice Phenomenology" 
Q^s already providing important contributions in particle 
l^hysics, and that it is developing into the premier quan- 
^^itative tool for non-perturbative field theory. For some 
QQphysical quantities lattice computations have been per- 
formed successfully for several years, and the emphasis 
^^s now on controlling and reducing the systematic errors 
^^in these computations. These quantities include the lep- 
^Sjtonic decay constants of mesons discussed in sections 2 
O^and 3 below. For other quantities, such as the evalua- 
' 'tion of the Isgur-Wise function of semi-leptonic decays 
Q^of heavy mesons discussed in section 4, lattice studies 
^^are only just beginning. 

Oh Before presenting some recent results, perhaps it is 
^^ecessary to explain briefly what lattice calculations 
QJare, and to discuss the major sources of uncertainty 
n these computations. The starting point for lattice 
studies is the evaluation of the functional integral: 
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where A and ip represent gluon and quark field respec- 
tively, 0{xi,X2, ■ ■ ■ ,Xn) is a multilocal operator com- 
posed of the fields, S represents the action and Z is the 
partition function. I{xi, X2, • • • , Xn) is the vacuum ex- 
pectation value of the operator O. The functional inte- 
gral I is evaluated by discretising space and time, and 
generating field configurations weighted by the Boltz- 
man factor e""^. The physical quantities which can be 
studied depend on the choice of operator O. For ex- 
ample, by choosing O to be a bilocal operator of the 
form 0{xi,X2) = Jh{x2)j\{xi), where j| and ,7^ are 
interpolating operators which can create or annihilate 
the hadron h, the propagation of hadron h is studied. 



allowing one to evaluate its mass and the matrix ele- 
ment {Q\Jh\h). By evaluating 3-point correlation func- 
tions with 0{xi,X2,xz) = Jh2{xz)X{x2)j\ {xi), where 
X is some local operator, we can evaluate matrix ele- 
ments of the form {h2\X{0)\hi) . For many fundamental 
quantities in particle physics, particularly in studies of 
hadronic structure and weak decay processes, the (non- 
perturbative) strong interaction effects can be expressed 
as operator matrix elements of this kind. Hence the im- 
portance of lattice simulations. 

I will attempt to suppress any discussion of the tech- 
nology of lattice computations, nevertheless it is diffi- 
cult for me to avoid using two pieces of lattice jargon. 
The first is /3 = 6/(jfo(a), where go is the bare cou- 
pling constant and a is the lattice spacing. It is conve- 
nient in lattice simulations to fix (3 (and hence go(a)), 
and to determine the corresponding value of the lattice 
spacing by comparing the lattice prediction for some 
physical quantity to its physical value, (rather than the 
apparently more natural procedure of fixing the lattice 
spacing and determining the corresponding value of go). 
The second piece of jargon is k, the Wilson hopping pa- 
rameter, which is proportional to the coefficient of ipip 
in Wilson's discretisation of the quark action (as well 
as generalisations of this like the SW action mentioned 
below). K is therefore a measure of the quark mass, 
and the "critical" value of k, which corresponds to zero 
renormalised quark mass (and zero pion mass), is called 

The numerical computation of the functional inte- 
gral in eq.(l) leads to an evaluation of operator matrix 
elements "from first principles". However, there are a 
number of approximations in these calculations, lead- 
ing to uncertainties in the final results. First of all we 
have the "statistical" errors, that is the errors due to 
the fact that we are estimating the functional integral 
by sampling the integrand at a finite number of field 
configurations. The size of these errors can be esti- 
mated, using standard statistical methods, by observ- 
ing how the result varies as configurations are added or 
removed. As will be clear from the results presented 



below, for many quantities (on typical lattices) it ap- 
pears that the statistical errors are adequately small 
for 50-100 field configurations (although there are some 
important exceptions to thisO. 

More problematical are the systematic errors. These 
include "finite volume effects" , i.e. errors due to the fact 
that the integrals are evaluated with space-time taken 
to be finite. These can be studied by repeating the 
simulations on lattices of different sizes, and this is now 
being done more frequently. Theoretically, for a range 
of interesting quantities, it is known that finite volume 
effects decrease exponentially with the volume [1], and 
numerically the effects appear to be small on currently 
used lattices, at least for quenched calculations. How- 
ever for this to be the case the numerical simulations 
are performed with light quark masses which are heav- 
ier than the physical ones (typically corresponding to 
pions with masses in the range of about 400 MeV - 
1 GeV), and the results are then extrapolated to the 
physical limit (which is essentially the chiral limit for 
which the quarks are massless). I mention in passing 
that the dependence of masses and energy levels on the 
spatial volume of the lattice can be used to measure 
scattering lengths [1]. 

A second source of systematic error is due to the 
finiteness of the lattice spacing. Again these errors can 
be studies by performing high statistics runs on lattices 
with different values of the lattice spacing (i.e. at differ- 
ent values of An example of this will be presented in 
sec. 2. Another approach is that of "improvement" [2] 
in which the discretisation errors are formally reduced 
through the use of an "improved" action and operators. 
Some of the results presented below have been obtained 
with the use of the fermion-action proposed by Sheik- 
holeslami and Wohlert (the SW or "clover" action) [3]. 
With this action, and with the use of improved opera- 
tors, the discretisation errors are reduced from those of 
0(a) (present in simulations with Wilson fermions) to 
ones of O(asa) [4]. An exciting possibility is that, by 
the use of the renormalisation group transformations, 
it may be possible to construct "perfect actions", for 
which the discretisation errors will effectively be elimi- 
nated, and yet which will be practicable for numerical 
simulations [5]. Whether this will be possible should 
become apparent during the next year or so. 

Finally there is the problem of "quenching" . The re- 
sults presented below have been obtained by neglecting 
the effects of quark loops (but including all gluonic ef- 
fects), and it is difficult in general to estimate the error 
this induces. Where the results of quenched calculations 
can be compared to experimental data, they are gener- 
ally in good agreement suggesting that these errors are 
moderate (an important exception may be the hyper- 
fine splittings in quarkonia). Sometimes it is possible to 



use the renormalisation group equations to relate two 
physical quantities, and the difference between results 
obtained by using zero fiavours or the physical number 
of fiavours in quark loops provides an estimate of the 
effects of quenching. However for the lattice computa- 
tions to be treated as a truly quantitative technique, the 
effects of quark loops should be included, for which bet- 
ter algorithms or theoretical developments (such as the 
perfect actions mentioned in the preceeding paragraph) 
are needed. Otherwise we will have to wait for several 
years (5-10 ?) while the improvements in computing 
resources become sufficiently powerful to cope with the 
evaluations of the fermion determinants, necessary for 
full QCD computations. 

In this short talk I clearly have to be selective in 
the subjects I cover. I hope nevertheless that the topics 
I have chosen will provide an accurate picture of the 
status of computations in lattice QCD, illustrating both 
what has been, and is being, achieved, and also some of 
the outstanding difficulties and attempts to overcome 
them. I will not have time to discuss many interesting 
lattice computations for field theories other than QCD. 

1 also regret that I will not be able to compare the lattice 
results with those from model calculations or QCD sum 
rules. 

2 Decay Constants of Light Mesons 

The leptonic decay constants of light mesons, (e.g. /,r 
and fp) are obtained from matrix elements which are 
among the simplest ones to compute using lattice QCD. 
They are defined by 

< 0|A^(0)|7r(p) >= Up, (2) 

and 

<0\V,(())\p(p)>=e,ml/fp (3) 

where is the polarisation vector of the p-meson, and 
A and V are the axial- vector and vector currents respec- 
tively. The experimental values of the decay constants 
are U = 132 MeV (U/mp ~ 0.17), and l/fp = 0.28(1). 

These decay constants have been computed for sev- 
eral years now. During the last year there have been a 
number of high statistics evaluations, and I will restrict 
my discussion to these studies. In particular the GFll 
Collaboration [6] have evaluated the decay constants at 
several values of the lattice spacing with sufficient pre- 
cision to be able to attempt an extrapolation to the 
continuum limit a = 0. In Fig.l I present the lattice 
results for j-^ j mpj Z a obtained by the CFll [6] and 
APE [7] collaborations, using Wilson fermions *. Za is 

*Tlie results presented liere for tlie GFll collaboration are 
ones whicli have been corrected since tlieir original preprint was 
circulated 



Figure 1: Values of j-^ j mpj Z a measured by the GFll 
and APE Collaborations. 



Group 


a 


a-^ GeV 


Lattice Size 


Configs 


GFll 


5.70 


1.36(2) 


16^ X 32 


219 


GFll 


5.93 


2.00(5) 


24^ X 36 


217 


APE 


6.00 


2.23(5) 


18^ X 32 


210 


GFll 


6.17 


2.78(5) 


30 X 32^ X 32 


219 



Table 1: Parameters of the Simulations by the GFll 
and APE collaborations whose results for the decay con- 
stants are presented in this talk 



the renormalisation constant relating the lattice axial 
current to the physical one, and is calculable in per- 
turbation theory. In addition to the directly measured 
values, I have also included in Fig.l, the results ob- 
tained by extrapolating these values to physical quark 
and pion masses. The results from the two collabo- 
rations are clearly in excellent agreement. The lattice 
sizes and numbers of gluon configurations used in each 
simulation are presented in table 1. The values of the in- 
verse lattice spacings presented in table 1 were obtained 
from the measured values of the mass of the p-meson 
from the corresponding simulation. 

The results in Fig.l show a surprisingly mild depen- 
dence on the lattice spacing. In order to extract the 
physical quantity fir /nip, they need to be multiplied by 
the renormalisation constant Za (which also depends 
mildly on the lattice spacing). Following the sugges- 
tion of Lepage and Mackenzie [8], I use the expres- 
sion Za ~ (1-0.31 6/(47r/3) (8kc)'*)/8kc, which includes 
the effects of a partial summation of tadpole diagrams 
(which are lattice artefacts giving large coefficients in 
lattice perturbation theory) [9, 10, 11]. The Z^'s ob- 





fw/mp 




5.70 


0.176(7) 


0.41(2) 


5.93 


0.165(5) 


0.36(2) 


6.00 


0.167(7) 




6.17 


0.165(5) 


0.34(2) 



Table 2: Values of fn/rrip and 1/ fp obtained by the 
GFll and APE collaborations. 

tained in this way are 0.67 at /3 = 5.7, 0.74 at /3 = 5.93, 
0.75 at P = 6.0 and finally 0.77 at /3= 6.17. The re- 
sults for fn/rrip obtained using these values of Za are 
presented in table 2. It must be stressed that the er- 
rors given in table 2 are statistical errors only. From 
these results we see that the dependence on the lattice 
spacing is remarkably small for this particular quantity, 
and that the results are in excellent agreement with the 
physical value. 

The situation is a little different for the decay con- 
stant of the /9-meson. Following the analogous proce- 
dure to that for f-n /nip above, we obtain (from the data 
of the GFll collaboration) the values of 1/ fp in the 
third column of table 2. If the behaviour with the lat- 
tice spacing is linear, the results for 1/ fp extrapolate 
to about l/fp = 0.25(3) at a = [6]. If this is the 
case, then clearly there are significant 0(a) effects in 
the evaluation of l/fp at around /3 ~6.0-6.2 for Wilson 
fermions, however further work is needed to establish 
that the behaviour with a is indeed linear all the way 
down to /3 = 5.7. 

Similar studies are beginning also with the SW- 
action, for which the discretisation errors are formally 
reduced [7, 12]. For example the UKQCD collabora- 
tion [12] find f^/mp = 0.14(1) at /3 = 6.2 and the APE 
collaboration find fTr/nip = 0.16(1) at /3 = 6.0. Over 
the coming months it will become possible to study 
the a-dependence in some detail. I should also men- 
tion that it is practicable, in simulations using the SW 
fermion action, to determine the renormalisation con- 
stants Za and Zy non-perturbatively, by imposing the 
chiral Ward Identities [13], removing one source of un- 
certainty. 

3 Decay Constants of Heavy Mesons 

In this section I will review the status of computa- 
tions of the decay constants of heavy-light pseudoscalar 
mesons, i.e. mesons with a heavy quark (anti-quark) 
and a light anti-quark (quark). The decay constant fs 
(in conjunction with the _B-parameter oi B°-B° mixing) 
is an unknown parameter needed for the determination 
of the elements of the CKM matrix, and the CP violat- 
ing phase in particular. I will present results obtained 



using two different approaches. The first method is an 
extension of the calculations described in section 2 but 
with the mass of the heavy quark in the region of that 
the charm quark (for quarks which are much heavier 
than this, the Compton wavelengths are smaller than 
the lattice spacing). I will refer to this method as sim- 
ulations with propagating heavy quarks. The second 
method involves simulations of the Heavy Quark Effec- 
tive Theory directly on the lattice, and I will refer to 
this method as simulations with static heavy quarks. I 
will not have time in this talk to describe lattice studies 
in heavy quark physics using the non-relativistic formu- 
lation of QCD (see however item c in section 5). 

Simulations with propagating heavy quarks have be- 
en performed for several years now, and for example, at 
the 1989 conference on Lattice Field Theory [14], Steve 
Sharp e summarised the results for fo as 

= 180±25±30MeV (4) 

Results from recent simulations all lie in this range. 

In the heavy quark effective theory the decay con- 
stant of a heavy pseudoscalar (P) meson is predicted to 
behave as a function of its mass as follows: 

fp = [(a^iMp j)-'^^" [1 + 0(a,)] + Oil/Mp)] 

_ (5) 

In lattice simulations it is the constant A which is eval- 
uated, following the method proposed by Eichten [15]. 
I will denote by f^^^ the value of fp obtained in this 
way, i.e. obtained by dropping the 0(1/Mb) terms. 
Early results for f^^^ gave surprisingly large values: 

^stat ^ 310±25±50MeV ref.[16] (6) 
/|tat ^ 366 ±22±55MeV ref.[17] (7) 

Clearly these results could only be made consistent with 
those in eq.(4) if there were large negative 0(l/Mp) 
corrections in the charm region, and significant ones for 
the 5-meson, so that fp ^ .fs"^^- It was therefore im- 
portant to check whether this was the case, and it was 
found in simulations with propagating quarks that in- 
deed the quantity fp/VMp{a,,{Mp) f/l^° does increase 
as Mp is increased [18, 19]. 

A compilation of the results for f^^^ is presented 
in table 3 [20]. In the last column I present the re- 
sults as they were quoted in the publications, but I also 
present the values of the lattice spacing and renormal- 
isation constant which were used to obtain the result. 
Part of the reason for the spread of results is due to 
different procedures, particularly in the choice of Z^*"*, 
but there is still some debate whether all of the col- 
laborations can isolate the contribution of the ground 
state sufficiently accurately. To make further progress. 



more work on improving the evaluation of Z^"" and 
above all on isolating the ground state is needed, and 
is in progress. There is some preliminary evidence that 
fstat j^^y |jg decreasing as a ^ [23, 24], and it will 
be very interesting to be able to check this when the 
precision of the calculations improves. 

This year two groups have presented new results ob- 
tained with propagating heavy quarks. Bernard, Lab- 
renz and Soni, have performed simulations at /3 = 6.3 
(20 configurations on a 24^ X 55 lattice) and at /3 = 6.0 
(19 configurations on a 16^ x 39 lattice and 8 config- 
urations on a 24^ X 39 lattice) using Wilson fermions. 
These authors attempt to reduce the systematic errors 
associated with O(mQa) effects (where Q is the heavy 
quark) by modifying the heavy quark propagator fol- 
lowing a procedure based on the structure of the free- 
field propagator [25, 26], (it will be very interesting to 
see whether this procedure for reducing O(mQa) effects 
in simulations with Wilson fermions can be established 
when quantum loops are included, and this work is in 
progress [27]). Bernard, Labrenz and Soni quote 

fp = 187 ± 10 ± 34 ± 15 MeV (8) 
fp = 208 ± 9 ± 35 ± 12 MeV (9) 

These values are based on their results with both prop- 
agating and static heavy quarks. 

The second collaboration to present results this year 
with propagating quarks is the UKQCD collaboration 
which performed simulations at /3 = 6.2 (60 configura- 
tions on a 24^ X 48 lattice), and /3 = 6.0 (36 config- 
urations on a 16^ X 48 lattice) using the SW fermion 
action. Since the leading discretisation errors for the 
SW action are of O(asmQa) (instead of O(mQa) as 
for Wilson fermions), it was important to check that 
the results found with Wilson fermions [18, 19], and 
the dependence on mq in particular, are reproduced ^ . 
In Fig. 2 I show the results for the "scaling" quantity 
fp/\fMp{as{Mp)fll^° as a function of l/Mp obtained 
by the UKQCD collaboration from their simulation at 
(3 = 6.2. The open points correspond to the measured 
values at three values of the mass of the light quark (de- 
creasing as the mass of the light quark is decreased) and 
the solid points are the results obtained after extrapo- 
lation to the chiral limit. The solid line is a linear fit 
to the heaviest three of the four points, and the broken 
line is a quadratic fit to all four points. The increase of 
this quantity as Mp inreases is manifest. Also on this 
figure is shown the value of /^^^^ obtained from 20 of 
the 60 configurations, which is about 2-3 standard devi- 
ations above the value obtained by extrapolation from 
the results with propagating quarks. From their results 

^ UKQCD have also analysed 18 (of tlie 36) configuratioiis on 
tlie lattice at /? = 6.0 using Wilson fermions, and compare tlie 
results for tlie two different fermion actions 



Ref. 


Action 


P 


GeV 


rjstat 
A 


/Stat jyjgv 


[21] 
[16] 
[17] 
[7] 
[7] 
[2i)\ 

[20] 
[22] 
[23] 


Wilson 
Wilson 
Wilson 
Wilson 
SW 

SW 
Wilson 
Wilson 


5.9 
6.0 
6.0 
6.0 
6.0 

6.U 

6.2 
6.3 

5.74, 6.0, 6.26 


1.75 
2.0 ±0.2 
2.2 ±0.2 
2. 11 ±.05 ±.10 
2.05 ±0.06 
2.0 _ 2 
2.7 1 I 
3. 21 ±.09 ±.17 
1.12, 1.88, 2.78 


0.79 

0.8 

0.8 

0.8 
0.89 

U. / 8 

0.79 
0.69 
0.71(8) 


319 ± 11 
310 ± 25 ± 50 
366 ±22 ±55 
350 ±40 ±30 

370 ± 40 

OSR +8+67 
/no _ 10 - 42 

9cq + 16 +105 

235 ±20 ±21 
230 ±22 ±26 



Table 3: Compilation of Lattice Results for /j 



I would like to conclude this section by briefly men- 
tioning a number of related quantities. First of all, as 
is clear from fig. 2 the decay constants decrease as the 
mass of the light quark decreases (this is interpreted 
as being due to the fact that the size of the meson in- 
creases, so that the wavefunction at the origin is re- 
duced). Lattice simulations typically give a result 10- 
20% larger for Jb, and Jb, than for Jb and Jb (the 
errors in the ratio are small), and for example UKQCD 
quote /£)^ =212 1 4 1 '^^MeV, to be compared to the 
experimental results of Jb, = 232 ± 45 ± 20 ± 48 MeV 
[29] and /b, = 344 ± 37 ± 53 ± 42 MeV [30]. 

Another important quantity is the (renormalisation 
group invariant) 5-parameter for B° — B° mixing for 
which the ELC collaboration find fB\/BB = 220 ± 
40 MeV and UI^BbM UIBb) = 1.19 ±0.10 



Figure 2: Values of fp/Ml/^{asiMp)/asiMB))'-'^^'^°^ 
measured by the UKQCD Collaboration. 

with propagating quarks UKQCD quote: 

fB = 160 1 ^ 1 II MeV (10) 

fB = 185 1 ^ + 42 MeV (11) 

The reason for such asymmetric errors in eqs.(lO) and 
(11) is the uncertainty in the value of the lattice spacing 
obtained from the physics of light hadrons, for which 
the UKQCD collaboration quote a"! = 2.7 1 [CeV. 

My summary of the lattice results for the decay con- 
stants of heavy mesons, based on the above and earlier 
simulations, are 

fB = 180 ±40 MeV (12) 
fB = 200 ± 30 MeV (13) 

An important step now will be to study the dependence 
of the results obtained with the SW fermion action on 
the lattice spacing, and this work is in progress [28]. 



4 The Isgur Wise Function 

Semi-leptonic decays of heavy mesons are an important 
set of processes for studies of the standard model, and 
in the determination of the elements of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix (particularly for the 
Vcb matix element). In the heavy quark effective theory, 
the two form-factors for the decay B ^ D + leptons 
and the four form-factors for decay B D* ± leptons, 
are all given in terms of a single unknown function of uj, 
(u) = vb-vb, where vb and vb are the four- velocities of 
the B and D or D* mesons) [31]. Moreover this func- 
tion, $(tt)), known as the Isgur- Wise function [31], is 
normalised at the zero-recoil point, ^(1) = 1. Experi- 
mental studies of these decays give results for iV^cil S.i'^) 
for values oi lc > 1, so a determination of Vcb requires 
an extrapolation of the experimental results to oj = 1. 
Lattice simulations can provide a determination of ^(oj) 
and also test whether the charm quark is sufficiently 
massive for the heavy quark effective theory to be use- 
ful for B ^ D and D* decays. 

A convenient way of determining ^(t<j) is by evaluat- 



Figure 3: A comparison of the ARGUS experimen- 
tal data for lycil'J(i^) with the lattice results from the 
UKQCD Collaboration. 

ing the elastic matrix element 

< Dip')\crc\Dip) >= (p + p'rf+iq^) (14) 

where q = p — p' and the single form-factor /"*", when 
considered as a function of uj (q^ = 2m\i{\ — uj) ), is 
equal to the Isgur-Wise function (once radiative correc- 
tions are included). The matrix element in eq.(14) can 
be evaluated using the techniques which have been de- 
veloped for decays of charmed mesons (i.e. for D K 
or TT and D K* or p semi-leptonic decays) [.32, 3,3]. 
Two groups have recently presented results for ^(t<j), 
Bernard, Shen and Soni [34] and the UKQCD collabo- 
ration [35, 36]. 

In fig. 3 I present the results from the UKQCD col- 
laboration, obtained at /3 = 6.2 from 60 configurations 
using the S'VF-action, and after extrapolation to the chi- 
ral limit for the light quark masses. Fitting the lattice 
results to Stech's relativistic-oscillator parametrisation 
[37]: 

^(.) = ^exp (-(2,^-1)^) (15) 

the UKQCD collaboration find p"^ = 1.2 t l- Also 
on Fig. 3 are the data points from the ARCUS exper- 
iment, and by comparing the lattice determination of 
i{uj) to the ARGUS data, the UKQCD collaboration 
find V"cJx/rB/1.48ps = 0.038 where the first 

set of errors is due to the experimental uncertainty and 
the second is due to the uncertainty in the lattice de- 
termination of p'^ . An extension of this work to another 
value of (3 is in progress. 



Bernard, Shen and Soni perform simulations with 
Wilson fermions, (but modifying the heavy quark prop- 
agator as discussed in section 3) at /3 = 6.0 (19 configu- 
rations on a 16^ X 39 lattice and 8 configs on a 24^ x 39 
one) and j3 = 6.3 (20 configurations on a 24^ x 61 lat- 
tice) They fix p' to be 0, and limit their calculations to 
the two lowest values of however from their simu- 
lations at two values of (3 and two values of the heavy 
quark mass they are able to obtain results for a range 
of w's. These authors quote p'^ = 1.41 ± 0.19 ± 0.19 
based on their results for the largest values of the light 
quark mass (an improvement in statistics is needed to 
control the extrapolation to the chiral limit). It would 
be very interesting to extend the range of momentum 
values in this work so as to be able to study better the 
systematic errors. 

So far S,(uj) has only been determined from the D 
D matrix element in eq.(14). A very interesting set of 
checks will be performed during the next few months to 
establish whether the expectations of the heavy quark 
effective theory for the relation between the form factors 
in Pseudoscalar Pseudoscalar and Pseudoscalar 
Vector transitions are satisfied, for heavy quark masses 
in the region of the mass of the charm quark. 

5 Some Other Lattice Studies in Particle Phys- 
ics Phenomenology 

I would like to end this talk by briefly mentioning a few 
other important studies in particle physics phenomenol- 
ogy which I haven't had time to discuss in detail. 

a) Determination of the Strong Coupling Constant: 
El-Khadra et al. have calculated the strong cou- 
pling constant by computing the 1S'-1_P mass split- 
ting in charmonium [38], finding 

a3j^(5 GeV) = 0.174 ±0.012 (16) 

Liischer et al., are attempting to compute a., by 
using the size of the lattice as the dimensionful 
parameter[39]. 

b) Potentials: The potentials between static heavy 
quarks can be computed accurately [40, 41] and 
used to determine a.,. 

c) Quarkonia: One of the exciting recent develop- 
ments has been the use of non-relativistic QCD for 
studies of the spectrum and properties of Quarko- 
nia [42, 43]. Kronfeld and Mackenzie are also try- 
ing to generalise this approach to construct a field 
theory which would be applicable for all quark ma- 
sses [25, 26], and it will be very interesting to see 
whether this can be achieved. 



d) Bk- The parameter Bk, which parametrises the 
strong interaction effects in K°-K° mixing has been 
caicuiated by severai groups during recent years. 
Lusignoii et ai. [44] summarise the iattice resuits 
as Bk = 0.8 ± 0.2, where Bk is the renormaii- 
sation group invariant 5if-parameter. The most 
precise resuits came from Kilcup et al. [45], using 
the staggered formulation of lattice fermions, who 
performed simulations at several values of /3 and 
extrapolated their results to zero lattice spacing. 
Recently (indeed since this talk was presented), 
Sharpe has demonstrated that the discretisation 
errors in these calculations are of O(a^) and not 
0(a) [46], considerably reducing the uncertainty in 
the extrapolation to the continuum limit. The au- 
thors of ref.[45] quote B^ = 0.825 ± 0.027 ± 0.023 
(preliminary). 

e) A7 =1/2 Rule: There has been no progress on 
this very fundamental problem recently. The prob- 
lem (for Wilson fermions) is how to subtract accu- 
rately the diveregence which occurs when the 
operators of dimension 6 which contain the strong 
interaction effects for this process, mix with the di- 
mension 3 operator sd (through the so called eye- 
diagram). Although the subtraction can be per- 
formed in principle, the final results have errors of 
0(100%). 

f) Semi-Leptonic Charm Decays: The semi-leptonic 
decays D K,t:,K* oi p -\- leptons have been 
studied for several years now [32, 33, 47], giving 
pleasing results both for the form-factors at zero 
momentum transfer, and for the behaviour of the 
form-factors with momentum transfer. This year 
the ELC Collaboration has presented results at 
/3 = 6.4 with Wilson fermions [48], and although 
these have larger statistical errors than some of the 
earlier ones, the results are consistent (e.g. for the 
D K decay ELC quote /+(0) = 0.65(18)). In 
ref.[48] use of the heavy quark effective theory is 
made to develop a method for extrapolating the 
results to those for 5-decays. 

g) & ^ 57: At this conference we have heard from 
the CLEO collaboration about their observation of 
the process B K*j [49]. In the standard model 
this process occurs through "penguin" diagrams, 
and the rate is sensitive to possible new physics. 
In order to determine whether the observed rate 
is consistent with the prediction of the standard 
model it is necessary to evaluate the strong in- 
teraction effects. This is being done by two lat- 
tice groups, Bernard, Hsieh and Soni [50] and the 
UKQCD collaboration [51], both using propagat- 
ing heavy quarks. In fig. 4 we see the measured re- 
sults from the TIKQCD collaboration for the form- 
factor Ti which contains the strong interaction ef- 



Figure 4: A comparison of the form-factor Ti deter- 
mined from the CLEO measurement and lattice com- 
putations. 

fects for this process. The open squares are results 
obtained assuming that Ti is independent of the 
mass of the light-quark, whereas the diamonds are 
the results obtained after extrapolation of the data 
to the chiral limit (there is no observed dependence 
on the light quark mass, but the errors grow sig- 
nificantly if an extrapolation is attempted). The 
crossed square and diamond are the correspond- 
ing points after linear extrapolation in the heavy 
quark mass to the physical value of niK* / Mp (Mp 
is the mass of the heavy pseudoscalar). The value 
deduced from the CLEO data, assuming that mj 
= 150 GeV (the dependence on nit is mild) is also 
shown, as is the extrapolated value from the simu- 
lation of Bernard, Hsieh and Soni (BHS). The lat- 
tice results suggest that the observed rate is con- 
sistent with the standard model contribution. 

6 Conclusions 

I hope that I have managed in this talk to demonstrate 
that lattice simulations are making important contribu- 
tions to particle physics phenomenology, and in particu- 
lar to many processes which are under intensive experi- 
mental investigation. I have also tried to explain a little 
about the systematic errors present in the calculations 
and about attempts to reduce and control them. Lattice 
simulations are developing into the major quantitative 
tool for non-perturbative strong interaction physics. 
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